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ABSTRACT 

By  using  the  Dugdale  model  for  a  crack  in  a  plate,  an  improved  formula  was  derived 
for  the  strain  energy  release  rate,  G.  The  formula  has  the  same  form  as  the  solution  for 
a  linear  elastic  plate,  except  a  correction  factor  is  used  which  corrects  for  both  the  effect 
of  yielding  and  the  finite  width  of  the  plate.  Curves  are  presented  giving  the  values  of  the 
correction  factor,  and  they  indicate  that  the  nominal  stress  to  yield  stress  ratio  has  a 
pronounced  effect  on  the  strain  energy  release  rate. 


111 


AFFDL-TR-65-186 


TABLE  OF  CONTENTS 

SECTION  PAGE 

I  INTRODUCTION  .  1 

II  PROBLEM  SOLUTION  .  2 

1.  Adoption  of  Dugdale’s  Model  .  2 

2.  Calculation  of  Axial  Rigidity  .  3 

3.  Strain  Energy  Release  Rate  .  5 

III  REFERENCE  .  7 

APPENDIX  I  Integrating  the  Crack  Boundary  Displacements  .  8 

APPENDIX  II  Calculating  the  Load,  P  .  10 


v 


AFFDL-TR-65-186 


ILLUSTRATIONS 

FIGURE  PAGE 

1  Model  of  Dugdale  Crack  .  22 

2  Variation  of  the  Strain  Energy  Release  Rate  Correction 

Factors  with  Ratio  of  Crack  Length  to  Strip  Width,  k  .  23 

3  Displacement  of  the  Crack  Boundary  .  24 

4  Dimensions  of  Centrally  Cracked  Plate  Under  Uniaxial 

Tension  .  24 


vi 


AFFDL-TR-65-186 


A 


A 


n 


B 

E 

G 

G 

K 


L 

P 

P' 

T 

U 

Y 


a 


c 

e 

e 

o 

k 

t 


v 


v 

a 


v 

o 


X,  y 


cr 


LIST  OF  SYMBOLS 

Crack  surface  area 

Plate  gross  cross  sectional  area 

Plate  net  cross  sectional  area 

Plate  width 

Young’s  modulus 

Strain  energy  release  rate 

Strain  energy  release  rate  at  onset  of  fast  fracture 

Plate  axial  rigidity.  Stress  intensity  factor 

Plate  length 

End  load  on  plate 

Total  load  across  net  section 

Applied  tension  stress  normal  to  crack 

Elastic  energy  in  plate 

Tensile  yield  stress 

Half  crack  length 

a  +  p 

Natural  base  of  logarithms 
Given  extension  of  plate. 

2a/B  ,  Ratio  of  crack  length  to  plate  width 
Plate  thickness 

Displacement  of  crack  boundary 
Displacement  at  crack  tip 

Displacement  at  center  of  crack 

Cartesian  coordinates 
Stress  in  y  direction 


vii 


AFFDL-TR-65-186 


LIST  OF  SYMBOLS  (Cont’d.) 
p  Yield  zone  length 

/3  7TT/2Y 

a  Arccosh  x  .  Strain  energy  release  rate  correction  factor 

c 

y  Strain  energy  release  rate  correction  factor 


viii 


AFFDL-TR-65-186 


SECTION  I 

INTRODUCTION 

In  fracture  mechanics  theory,  G  is  defined  as  the  plane  stress  value  for  the  strain  energy 

release  rate  with  crack  extension  per  unit  of  crack  surface  area.  At  onset  of  unstable 

crack  propagation,  the  critical  strain  energy  release  rate,  G  ,  is  approximately  equal 

c 

to  the  plastic  energy  dissipation  rate  in  the  yielded  region  at  the  crack  tip.  The  size  of  the 
yielded  region  is  assumed  to  be  so  small  that  the  elastic  stress  distribution  is  unaltered 
by  the  presence  of  the  yield  zone.  Then  the  strain  energy  release  rate  can  be  calculated 
by  linear  elastic  analysis  methods.  Futhermore,  Gc  is  assumed  to  be  a  material  parameter 
governing  the  fracture  toughness  of  the  material  or  the  resistance  to  crack  extension. 


Unfortunately,  in  actual  practice,  G  has  been  found  to  vary  with  specimen  dimensions 

c 

and  with  nominal  (gross)  stress  levels.  Lorenz  (Reference  1)  found  a  definite  effect  of 
nominal  test  stress  levels  upon  measured  fracture  toughness  of  centrally  cracked  panels. 
Kraft,  Sullivan,  and  Boyle  (Reference  2)  introduced  the  concept  of  crack  extension  re¬ 
sistance,  R,  and  discussed  the  variation  in  R  with  specimen  width  and  crack  length. 

To  progress  with  the  use  of  G  in  fracture  mechanics  problems  it  is  necessary  to  derive 
a  more  accurate  solution  for  the  strain  energy  release  rate  which  takes  into  account  the 
plastic  deformation  at  the  crack  tip.  The  current  practice  to  correct  for  the  effect  of  yield¬ 
ing  is  to  arbitrarily  assume  an  increased  crack  length  which  extends  into  the  yield  zone. 
The  yield  zone  size  is  determined  by  applying  a  yield  criteria  to  the  linear  elastic  stress 
field  equations.  The  object  of  this  report  is  to  derive  a  more  reasonable  solution  for  the 
strain  energy  release  rate  by  using  a  model  which  does  not  assume  that  the  elastic  stress 
distribution  is  unaltered  by  the  presence  of  the  yield  zone.  This  is  accomplished  by  using 
the  Dugdale  (Reference  3)  model  for  which  the  yield  zone  size  is  formulated  from  equilibrium 
considerations.  With  the  Dugdale  model,  G  is  obtained  by  using  a  method  derived  by  Irwin 
and  Kies  (Reference  4).  A  tension  load,  P,  is  assumed  to  be  applied  to  the  plate,  and  a 
strain  energy  function  is  derived  from  Greenspan’s  (Reference  5)  formulation  for  the  axial 
rigidity  of  the  plate.  The  strain  energy  function  is  differentiated  with  respect  to  the  crack 
surface  area  to  obtain  the  strain  energy  release  rate. 
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SECTION  II 


PROBLEM  SOLUTION 


1.  ADOPTION  OF  DUGDALE’s  MODEL 


Dugdale’s  model  for  a  crack  in  tension  under  plane  stress  is  shown  in  Figure  1(a).  The 
model,  derived  from  the  mathematical  developments  of  Muskhelishveli  (Reference  6),  is 
based  on  the  following  assumptions: 

1.  Yielding  occurs  in  a  narrow  wedge  shaped  zone. 

2.  The  material  in  the  zone  is  under  a  uniform  tensile  yield  stress,  Y. 

3.  A  Tresca  yield  criterion  is  obeyed. 

4.  The  material  outside  the  zone  is  elastic  and  bounded  internally  by  a  flattened  ellipse 
of  length  2(a  +p )  where  a  is  the  half-length  of  the  crack,  and  p  is  the  length  of  the  plastic 
extension. 

5.  The  length  p  is  such  that  there  is  no  stress  singularity  at  the  ends  of  the  flattened 
ellipse. 

By  determining  the  plastic  zone  size  in  equilibrium  with  the  applied  stress,  Dugdale 
obtained  the  following  solution  for  the  plastic  extension: 


p  =  a  ( sac  0  -| ) 


(  I) 


where  /3  =  «jp  —  and  T  is  the  applied  tension  stress  normal  to  the  crack. 

Agreement  between  the  observed  and  calculated  values  of  p  has  been  exceedingly  good. 
Dugdale’s  experimental  results  showed  good  agreement  with  theory  Goodier  and  Field 
(Reference  7)  referred  to  numerous  experimenters  who  have  observed  the  Dugdale  type 
yield  zone.  Hahn  (Reference  8)  has  obtained  good  agreement  with  theoretical  p  in  tests 
on  thin  plates  of  silicon  steel  in  which  the  yield  zones  appeared  as  those  shown  in  Figure 
1(b).  Hahn  (Reference  9)  also  determined  that  the  yield  zone  configuration  begins  to  ap¬ 
proach  a  narrow,  tapered  Dugdale  type  when 


(2) 


p  ~  4t 


where  t  is  the  plate  thickness.  Finally,  this  writer  has  observed  close  agreement  with 
Equation  1  in  tests  of  0.02  inch  thick  sheets  of  AM350  and  AM355  steel  in  which  the  yield 
zones  were  also  similar  in  shape  to  the  one  shown  in  Figure  1(b). 

Hahn  (Reference  9)  programmed  Dugdale’s  stress-field  solution  for  a  computer  and 
determined  that  some  terms  in  the  solution  were  negligibly  small.  Ignoring  these  terms, 
the  stress  gradient  is  described  by  the  equation 


arc  tan 


(3) 


2 


AFFDL-TR-65-186 


where  cr  is  the  stress  in  the  y  direction, 

a  -  arc  cosh  — —  and  c  =  a  +  p 
c  r  • 

The  displacement  of  the  crack  boundaries  has  been  derived  by  Goodier  and  Field  (Ref¬ 
erence  7)  for  the  Dugdale  model  and  is  given  by  the  equation 


*'  7T  [cos  ( 


sin  ( yS- 8  ) 
sin2  (  0  +  6  ) 


j  +  cos  $  -In  ^ 


(sinyS  4*  sin  G  ) 
(sinyS  —  sin  0  ) 2 


)] 


(4) 


where 

Poisson’s  ratio  is  taken  as  1/3. 
v  is  the  displacement  in  the  y  direction. 
E  is  Young’s  modulus. 

8  -  arc  cos  — 


At  the  center  of  the  crack  (x  =  O),  Equation  4  reduces  to  the  form 

T  o  n  /  sln^  +  l  \2 

*•  *  T  TF  M  ) 

At  the  tip  of  the  crack  (x  =  a)  ,  Equation  4  reduces  to 


v 

a 


E 


2  a 
0 


Jn 


sec/8 


(5) 


(6) 


Hahn  (Reference  9)  experimentally  measured  v  ,  and  found  good  agreement  with  Equation 

a 

6. 

Equations  3,  4,  5,  and  6  make  up  the  solutions  for  the  stress  gradient  and  displacements 
which  are  required  in  the  next  section  to  derive  the  axial  rigidity  expression. 

2.  CALCULATION  OF  THE  AXIAL  RIGIDITY 

Using  the  reciprocal  theorem  in  elasticity  on  a  perforated  and  an  unperforated  plate, 
Greenspan  derived  the  following  expression  for  the  axial  rigidity  of  a  perforated  plate 
subjected  to  uniaxial  tension: 


I 

K  =  - 1 -  (7 

I  +  —  E—  J  Up  dx 

PL  J  p 

where 

the  coordinates  x  and  y  are  in  the  width  and  length  directions,  respectively,  as  shown  in 
Figure  1(a). 
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K  is  the  axial  rigidity  defined  as  the  ratio  of  the  overall  extension  of  the  imperforated 
plate  to  that  of  the  perforated  plate, 

P  is  the  uniaxial  tensile  force  applied  to  the  plate. 

L  is  the  length  of  the  plate. 

u  is  the  y-direction  displacement  along  the  perforation  which  is  integrated  along  the 
P 

hole  boundary. 


As  shown  in  Equation  7,  and  load,  P,  and  the  integral  of  the  displacements  along  the  per¬ 
foration  must  be  known  to  determine  the  axial  rigidity,  K.  Using  Equation  3  to  calculate 
P  and  Equations  4,  5  and  6  to  calculate  the  integral  of  the  displacements,  the  following 
expression  for  the  axial  rigidity  was  obtained:  * 


K  = 


l  + 


C  (v  ) 


2  L 


nz  C(n ) 


where 


Cl.).  — r  JU  ..eg  4-  (  •”  •“  p  )  A  1 

TL  '  4p  'sin/9—  I'-* 


■[-(■ 


4  jUi 


m 


Z  „  I 


A  /  sin/34-1  ^2 
^  '  sin/9  -I  > 


)  ] 


72 


(8) 


N  =  n  +  S —  m z 


m  *  sec  /9 


B  . 


n  =  ^  is  the  ratio  of  plate  width  to  total  crack  length. 


(a) 


For  0  <  —  <  0.5, 


C  (  n  )  = 


(m-l  ) 


Y  (n  -m )  m  tan  2/9 

-  -  4-  -  +  - - — 

n  T  n  2n  $ 


sin‘  $ 


m 

~N~ 


Vcos  20 


JU 


( N  —  m  -J co  s  2  j9  )  (14-  i/co  s  2  /9  ) 
(N  +  m  ~y~coi  2  j9  )  ( I  —  V cos  2/9  ) 


(b  )  For  0.5  <  —  <  I, 


C  (n  )  = 


(m-l)  Y 


(n-m  )  m  tan  2  /9 


2n/9 


L  N 


2  sin2;9 
a/- cos2  0 


arc  tan 


(N-m  )  cos  20 
N-m  cos  2/3 


)] 


*  Some  details  of  this  computation  are  shown  in  Appendixes  I  and  II. 
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3.  STRAIN  ENERGY  RELEASE  RATE 

If  a  given  tension  load  is  applied  to  a  centrally  cracked  plate,  the  stored  elastic  energy 
is  inversely  proportional  to  the  axial  rigidity,  K.  Following  a  development  similar  to  Irwin’s, 
it  is  necessary  to  assume  an  initial  load  Pq  which  would  prevail  for  a  given  extension 

e^  of  the  plate  if  the  central  crack  did  not  exist.  Thus,  for  a  plate  of  thickness  t, 

Pn  L  PL 

e  = - -  -  ( 9) 

0  EBt  KEBt 


where  P  is  any  tensile  load  supported  by  a  plate  containing  a  transverse  central  crack. 
For  plane  stress,  the  elastic  energy  is 

1  P2L  I  2  EBt 

U  =  — - -  ' 

2  KEBt 

Letting  k  =  1/n  =  2a/B,  then  for  fixed  ends, 

dU  _  j_ _ 

dk  2  8°  L  dk 

Or,  since  the  crack  surface  area,  A,  is  equal  to  Btk, 

dU  I  dU  _  J_  T2  L  dK 
dA  Bt  dk  "  2  E  dk 


2  o  (_ 


2  EBt  dK 


where  T  =  P  /  (Bt). 

Finally,  using  the  expression  for  K  in  Equation  8,  the  strain  energy  release  rate  for  the 
Dugdale  model  is 


7T  T  2  a 

C(  v ) 

r  tan 

[l+  Cln,+  /JLN 

*  I 
-i)  J 

h 

E 

2  1 

rc”l'-cu’j 

r 

If  L»B  which  is  a  condition  satisfied  in  many  instances.  Equation  13  simplifies  to 


where 


dU 

-  =  G 

dA 


2 

IT  T  o  2 

~ —  r 


r 


C  (  v  ) 
2C(n)? 


[ 


I  +  C (n  )  + 


ton  /3  I 
/3(nN  -I  )  J 


(14) 


Equation  14  has  the  same  form  as  the  solution  for  a  linear  elastic  plate,  except  that 
y  is  a  correction  factor  which  corrects  for  both  the  finite  width  of  the  plate  and  the  effect 
of  the  plastic  deformation  at  the  crack  tip.  To  aid  in  using  the  correction  factor,  values 
of  y  are  listed  in  Table  1  for  appropriate  ranges  of  T/Y  and  k. 


To  show  the  significance  of  Equation  14,  the  equation  can  be  compared  with  the  extended 
crack  solution  reported  in  Reference  10.  The  extended  crack  solution  can  be  written  as  follows: 
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K 

E 


2 


G 


(15) 


where  in  this  equation,  K  is  the  stress  intensity  factor  expressed  as 


K 


2 


2  2  2 

For  a  first  approximation,  if  it  is  assumed  that  K  =ttT  aa,  Equation  15  becomes 


(16) 


G 


-  a 

E 


where 


a 


2  2 


tan 


{ 


7T  k 


2 


[ 


m!  ] } 


For  comparison,  the  values  of  a  and  y  are  plotted  in  Figure  2  for  several  ratios  of  T/Y.  This 
figure  shows  that  at  higher  ratios  of  T/Y,  the  extended  crack  solution  gives  significantly 
lower  values  for  G  than  the  improved  solution  based  on  the  Dugdale  model.  Also  of  particular 
interest  is  the  difference  in  a  and/ when  T/Y  is  zero.  The  reason  for  this  is  because  a  is 
a  solution  for  equally  spaced  collinear  cracks  in  an  infinite  plate.  If  the  values  of  y  for 
T/Y  equal  to  zero  are  compared  with  Dixon’s  (Reference  11)  correction  factor  for  a  finite 
width  plate,  the  agreement  is  excellent. 
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APPENDIX  I 

INTEGRATING  THE  CRACK  BOUNDARY  DISPLACEMENTS 

Due  to  complications  in  integrating  Equation  4  directly,  an  indirect  but  accurate  method 

of  integration  was  used.  If  v  is  plotted  along  the  axis  of  the  crack  from  x  =  -a  to  x  =  a, 

the  resulting  curve  has  an  elliptical  shape.  This  curve  can  be  generally  approximated 

within  three  significant  figures  by  the  following  equation  in  which  f  is  determined  from  the 

condition  that  v  =  v  at  x  =  a. 

a 


o 


where 


vq  is  obtained  from  Equation  5. 


v&  is  obtained  from  Equation  6. 


f  = 


l  - 


4  Qm  sec  /3 


'll 


JU  (-‘I? %+.!)* 

'  sin/S  -  i 

The  area  of  the  shaded  segment  of  the  ellipse  shown  in  Figure  3  gives  the  value  of  the 
desired  definite  integral.  The  integral  is 


t 

; 


-a 


vdx  =  av„  -I-  f arc  sin — 
a  o  f 


2  T 

7 T  a  -  C  (  v  ) 

E 


(18) 


where 


c ( w )  ...jb  +  rv"1  ”■■).£.(  *‘"f -1  j  ] 

T  L  \  4p  /  \  sln/9  -  I  /  J 


4  Jin  sec  (3  \* 

>  (1 Wyf 


Equation  18  is  the  solution  for  an  infinite  plate  with  the  correction  factor  C(v)  accounting 
for  the  affect  of  the  yield  deformation.  For  comparison,  values  of  C(v)  are  listed  below 
for  values  of  T/Y  ranging  from  0  to  0.8. 

T/Y  C(v) 

0  1 


0.2 

1.028 

0.4 

1.115 

0.6 

1.288 

0.8 

1.661 
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As  just  mentioned.  Equation  18  is  for  a  crack  in  an  infinite  plate  and  then  should  probably 
not  be  used  for  k  >  0.5. 


10 


AFFDL-TR-65-186 


APPENDIX  II 

CALCULATING  THE  LOAD,  P 

Using  a  method  similar  to  Greenspan’s,  the  load,  P,  can  be  determined  by  integrating 
the  axial  stresses  along  the  x  axis  from  x  =  0  to  the  plate  edge,  x  =  na.  The  symbols  used 
for  the  derivation  are  shown  in  Figure  4.  The  total  load  carried  by  the  plate  is 


where 


no 

p'  =  2f  J  a y  dx 


=  2Y  t  a  ( m  —  I  )  +  2 1  /  cry  dx 
ma 


/ 


T  ,  *in  2)3  ^ 

cr  =  T  +-  -  arc  tan  [  — — -  ] 

y  a  \  za  p  I 

H  €  “COS  2  p 


(19) 


2a 

6 


=  nr  a 


A  useful  simplification  in  integrating  the  stresses  can  be  obtained  by  making  the  following 
assumption: 


na 


J  arc  tan/ - 

ma  ,  ^  (  x  + 


c  s 


in  2J3 


dx 


-/x2  -  c2  cos  2)8  )  2  —  c* 


COS 


2/9 


(20) 


na 


c *  3 i  n*  j9  d  x 

( x  +  -v/x2  —  c*  cos  2)8  )*  -  c2  cos2)8 


This  simplification  was  checked  numerically,  and  for  the  ranges  required  for  T/Y  and  k, 
the  maximum  error  that  resulted  was  less  than  3  percent. 

The  integration  can  be  performed  by  making  the  following  change  in  variables: 

w  =  x  V  x 2  -  c 2  ,  dx  =  (  I  —  ~T~)  dw  (21) 

Since  cos2 /3  is  positive  for  O  <  T/Y<0.5  and  negative  for  0.5<T/Y<1,  two  solutions  are 
required  for  the  integration.  After  performing  the  integration,  the  total  load  is  as  follows: 

(  a  )  For  0  <  y-  <  0.5  , 


P  =  2  at  T 


(-^-)  (m-l)f  2  atT(  n-m) 
2otmT  tan  2)3  f  m  sin*)3 


2/9 


L  N 


Vcos  2)8 


Jk( 


(N-m  -/cos  2/8  )  ( I  +  Scot  Z  (3  ) 
(N  + m  )  ( I  - ) 


)] 


(22) 
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(b)  For  0.5  <  — —  < 


P*  =  2atT(“— j  (m-l  )  +■  2atT  (n-m) 


2atmT  ta n2)3 
2)3 


L  N 


2  sinZ 

^“C03  2)9 


^  arc  tan 


m  V“  cos  2/3 


—  arc  tan 


V-C03  2/3 


)] 


(23) 


where 


N  -  n  +-  -m? 


But,  the  end  load,  P,  may  be  substituted  for  the  total  load  across  the  net  cross  section, 
p'  ,  if 


A  g  A„ 


2ot 


(24) 


where 

is  the  plate  gross  cross  sectional  area  at  y  =  L/2. 
Ar  is  the  plate  net  cross  sectional  area  at  y  =  O. 
Making  this  substitution,  the  total  load,  P,  is 


P  =  TA  n  C  (n  ) 


(25) 


where 


(a)  For  O  <  <  0.5, 


C(n)  =  inili 1  1  +  (n-m)  +  m  tan  2  $  F(  _  _  ^in2  /3  (N-mVcos  2/3  )  (  1  + Vc os  2)3  )  j 

n  ^  n  2n£  L  ^  VcoTi^  (NfmVcos2yS  )  ( I-  \/cos  2/S ) 


(b)  For  0.5<  —  <  I  , 

y 


x  (m-l)  Y  .  (n-m)  ,  m  tan  2  /8  f.  m  2sin*i8  ( N  -  m)  V~cos  2/3  ~| 

C(n  )  = - —  -I - +-  — — 7}L —  1 - - - ,  =  ..CI.t:  arc  tan  - 

n  T  "  2n£  L  N  V-  co,  2)9  N-m  J 


N-m  cos  2  £ 
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When  T/Y  =  O,  the  expression  for  C(n)  reduces  to  the  form: 


C  (n  ) 


(26) 


It  is  interesting  to  compare  Equation  26  with  Greenspan’s  solution  for  C(n)  used  in  Refer¬ 
ence  4  to  obtain  the  strain  energy  release  rate.  This  solution  is 


c  (n  )  =  i - k* 


(27) 


For  equal  values  of  k.  Equation  27  gives  slightly  lower  values  for  C(n)  than  does  Equation 
26.  The  reason  for  this  is  that  Equation  27  was  actually  derived  for  a  circular  hole  in  a  plate 
and  Equation  26  was  derived  for  an  elliptical  hole. 
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Figure  1.  Model  of  Dugdale  Crack 

(a)  The  Dugdale  Model 

(b)  The  Actual  Crack  (Hahn’s  Experiment  Reference  8  ) 
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STRAIN  ENERGY  RELEASE  RATE  CORRECTION  FACTOR,  y  AND  a 


AFFDL-TR-65-186 


Figure  2.  Variation  of  the  Strain  Energy  Release  Rate  Correction 
Factors  With  Ratio  of  Crack  Length  to  Strip  Width,  k 
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Figure  3.  Displacement  of  The  Crack  Boundary 
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Figure  4.  Dimensions  of  Centrally  Cracked 
Plate  Under  Uniaxial  Tension 
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